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1.  Publications 


The  following  scientific  publications  have  resulted  from  the  research  supported  by 
Grant  AFOS R-9  1-026 1 : 

Competing  sources  of  heavy  C02  enhancement  in  the  stratosphere  and  magnitude 
of  heavy  03  enhancement 
Planet.  Space  Sci.  39,  945  (1991) 
by  D.  R.  Bates. 

Rate  coefficients  for  N*(u)  dissociative  recombination 
Planet.  Space  Sci.  39.  1297  (1991) 
by  D.  R.  Bates  and  J.  B.  A.  Mitchell. 

Dissociative  recombination  of  polyatomic  ions 

J.  Phys.  B:  At.  Mol.  Opt.  Phys.  24,3267  (1991) 
by  D.  R.  Bates. 

Single- electron  transitions  and  cluster  ion  super  dissociative  recombination 
J.  Phys.  B:  At.  Mol.  Opt.  Phys.  25,  3067  (1992) 
by  D.  R.  Bates. 

Association  of  CH^"  and  tt2:  conflicting  results 
Astrophys.  J.  375,  833  (1992) 
by  D.  R.  Bates. 

Nightglow  emissions  from  oxygen  in  the  lower  thermosphere 
Planet.  Space  Sci.  40,  211  (1992) 
by  D.  R.  Bates. 
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Emission  of  the  forbidden  red  and  green  lines  of  atomic  oxygen  from  the  nocturnal 
F-region 

Planet.  Space  Sci.  40,  893  (1992) 
by  D.  R.  Bates. 

Photoionization  of  excited  states  of  atomic  nitrogen 

J.  Phys.  B:  At.  Mol.  Opt.  Phys.  25,  1209  (1992) 
by  K.  L.  Bell.  K.  A  Berrington  and  C.  A.  Ramsbottom. 

R-matrix  calculations  of  low-energy  electron  scattering  by  oxygen  molecules 
Phys.  Rev.  Lett.  68,2011  (1992) 
by  C.  J.  Noble  and  P.  G.  Burke. 

Dissociative  recombination  when  potential  energy  curves  do  not  cross 
J.  Phys.  -B:  At.  Mol.  Opt.  Phys.  25,  5479  (1992) 
by  D.  R.  Bates. 

Photoionization  cross  sections  for  atomic  oxygen 
Planet.  Space  Sci.  40,  1419  (1992) 
by  K.  L.  Bell  and  R.  P.  Stafford. 

Stratospheric  heavy  ozone:  the  symmetric  isomer 
Planet.  Space  Sci.  40,  1573  (1992) 
by  W.  L.  Morgan  and  D.  R.  Bates. 
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Tidal  termolecular  tonic  recombination 

J.  Phys.  B:  At.  Mol.  Opt.  Phys.  25,  5421  (1992) 
by  VV.  L.  Morgan  and  D.  R.  Bates. 

Enigma  of  Yif  dissociative  recombination 
Planet.  Space  Sci.  41,  9  (1993) 
by  D.  R.  Bates,  M.  F.  Guest  and  R.  A.  Kendall. 

Dissociative  recombination:  Unexpected  twist  in  a  long  tale 

Comments  on  Atomic  and  Molecular  Physics  29,  53  (1993) 
by  D.  R.  Bates. 

Vibrational  excitation  of  products  of  dissociative  recombination 
Mon.  Not.  R.  Astron.  Soc.  263,  369  (1993) 
by  D.  R.  Bates. 

Cause  of  terrestrial  nightglow  continuum 

Proc.  R.  Soc.  Lond.  A  443,  227  (1993) 
by  D.  R.  Bates. 

Prevalence  of  rapid  dissociative  recombination  in  absence  of  crossing  of  potentials 
Proc.  R.  Soc.  Lond.  A  443,  257  (1993) 
by  D.  R.  Bates. 
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The  ionosphere 

In  “Relating  Geophysical  Structures  and  Processes’,  (American  Geophysical 
Union,  1993),  pp  11  -18. 
by  D.  R.  Bates. 

Exchange  processes  in  the  scattering  of  polarized  electrons  by  oxygen  molecules 
J.  Phys.  B.  At.  Mol.  Opt.  Phys.  27,  185  (1994) 
by  C.  Fullerton,  G.  Woste,  D.  G.  Thompson  and  K  Blum. 

Low-energy  electron  scattering  by  oxygen  molecules 

J.  Phys.  B:  At.  Mol.  Opt.  Phys.  (in  press)  (1994) 
by  K.  B.  Higgins,  C.  J.  Noble  and  P.  G.  Burke. 

Ionization  and  simultaneous  excitation  of  helium  atoms  by  electron  impact 
J.  Phys.  B:  At.  Mol.  Opt.  Phys.  (in  press)  (1994) 
by  A.  Raeker,  K.  Bartschat  and  R.  H.  G.  Reid. 
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2.  Dissociative  recombination 

Laboratory  afterglow  measurements  of  Zipf  (1980)  that  have  been  presumed  to  show 
that  the  recombination  coefficient  of  N^"(u)  is  independent  of  the  vibrational  level  v 
have  been  analysed  (Bates  and  Mitchell,  1991)  to  make  allowance  for  the  coupling  due 
to  symmetrical  resonance  charge  transfer 

N  +  (u)  +  N2(0)  ^  N  +  (0)  +  N2(y)  (1) 

and  also  for  the  vibrational  deactivation  in  collisions  with  Ne  atoms  present  in  the 
afterglow.  An  experimental  value  for  this  vibrational  deactivation  coefficient  that  has 
been  reported  (Dobler  et  a 1  ,  1982)  is  much  higher  than  would  be  expected  or  than 
can  be  reconciled  with  the  data  of  Zipf.  It  is  therefore  rejected.  The  analysis  shows 
that  the  recombination  coefficient  for  N^(0)  is  about  2.6  X  10_7cm3s~ 1  at  300  K  and 
that  the  recombination  coefficients  for  N^"(l)  and  N^(2)  are  much  smaller.  These 
coefficients  conflict  with  the  dissociative  recombination  cross  sections  versus  energy 
curve  obtained  by  the  merged  beam  method  (Noren  et  a 1  ,  1989).  Reasons  were 
presented  for  believing  that  the  curve  cannot  be  correct. 

Experimental  data  on  products  of  the  dissociative  recombination  of  polyatomic 
ions  have  been  discussed  and  rationalized  in  terms  of  the  valence  bonds  concerned 
(Bates.  1991b).  By  combining  the  data  with  theoretical  considerations,  branching 
ratios  for  the  various  dissociative  channels  have  been  derived.  There  are  two  main 
classes  of  saturated  polyatomic  ions:  those  for  which  the  number  of  valences  of 
the  ionized  atom  is  decreased  by  neutralization,  and  those  for  which  the  number  is 
increased.  The  hydride  ions  H30  +  and  H4N  +  are  simple  examples  of  the  valence 
decrease  class.  We  have  reasoned  that,  while  H30  +  may  shed  two  H  atoms  in 
dissociative  recombination  with  thermal  electrons,  H4N  +  cannot  do  likewise  because  of 
a  difference  in  the  OH  and  NH2  relic  radicals.  We  predict  that  a  bond  between  a  pair  of 
neutral  atoms  in  an  ion  like  CH30"*"H2  cannot  be  broken  by  dissociative  recombination. 
Although  the  number  of  valences  of  an  isolated  C+  ion  increases  on  neutralization,  we 
predict  that  the  dissociative  recombination  of  polyatomic  ions  containing  C+  normally 
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takes  place  through  a  radiationless  transition  that  causes  a  decrease  For  example 
when  CH^  sheds  an  II  atom  by  dissociative  recombination  the  CH.>  relic  radical  should 
be  mainly  in  the  *Aj  metastable  state  with  the  C  atom  having  two  paired  orbitals 
(  a  pure  “lone  pair")  so  that  it  is  effectively  divalent.  In  contrast  to  CH30+H.> 
dissociative  recombination,  we  expect  that  CH3C'+H2  dissociative  recombination  may 
break  a  bond  between  neutral  atoms.  In  this  event  ethylene  would  be  formed  From  an 
examination  of  the  distribution  of  the  measured  dissociative  recombination  coefficients 
of  polyatomic  ions  we  have  inferred  that  they  are  constrained  by  having  an  upper  limit, 
and  we  have  shown  that  this  may  be  understood  by  application  of  the  maximum  cross 
section  theorem  of  collision  theory. 

Energy  considerations  show  that  dissociative  recombination  of  the  proton-bridge 
ions  H20  •  H+  •  H20  and  N H 3  •  H +  •  N H  3  can  occur  only  through  single-electron 
transitions.  The  nature  of  a  single-electron  transition  that  leads  to  dissociative 
recombination  has  been  discussed.  Calculations  have  demonstrated  that  the  rate 
of  the  radiationless  transition  involved  is  extremely  fast  (Bates.  1992a)  We  have 
also  discussed  the  recombination  of  the  related  cluster  ions  H30+(H20)n  and 
H4N+(NH3)n  and  explained  the  experimental  results. 

Direct  and  indirect  dissociative  recombination  have  been  compared  If  the 
attractive  potential-energy  curve  of  the  ground  state  of  the  ion  is  not  crossed  bv  a 
suitable  repulsive  potential-energy  curve  of  the  molecule  formed  by  neutralization, 
the  indirect  process  can  be  much  faster  than  the  direct  process.  However,  in  the 
case  of  HCO+  (for  example),  its  rate-coefficient  a(I)  is  several  orders  of  magnitude 
less  than  the  measured  value  even  though  the  radiationless  transition  involved  is 
of  the  extremely  rapid  single-electron  type.  A  new  process,  multi-step  indirect 
dissociative  recombination,  that  entails  the  system  passing  successively  through  several 
intermediate  states  in  Ai/  =  1  steps,  has  been  introduced.  Examination  of  the  formula 
for  its  rate-coefficient  Q(I*),  obtained  from  high-order  perturbation  theory,  shows  that 
the  final  pre-d issociative  step  is  rate-limiting.  Extensive  computations  have  been  done 
on  the  overlap  integrals  between  the  bound  vibrational  wavefunctions  of  the  ion  and 
the  wavefu n c tion  for  the  separating  fragments  of  the  molecule  It  lias  been  concluded 
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from  these  that  multi-step  indirect  dissociative  recombination  provides  a  solution  to 
the  enigma  posed  by  the  laboratory  data  on  HCO  +  recombination  The  multi-step 
mechanism  is  the  key  to  understanding  dissociative  recombination  in  cases  where  the 
relevant  potentials  do  not  cross  and  the  transition  from  the  continuum  is  single-electron 
(Bates  I992e). 

The  high  dissociative  recombination  coefficient  measured  for  Hj  in  its  ground 
vibrational  level  has  appeared  irreconcilable  with  theory,  there  being  no  crossing  of 
potential-energy  curves  favourable  to  dissociative  recombination  at  low  temperature, 
h  has  been  shown  that  the  new  mult-siep  process  can  explain  the  experimental 
results.  Strong  support  for  the  mr'M-step  process  has  been  provided  by  measurements 
on  the  recombination  coefficient  which  is  comparable  in  magnitude  with  the 

recombination  coefficient  despite  the  quanta!  tunnelling  through  a  subslan.ial 
potential  barrier  that  is  involved.  In  earlier  theory,  such  tunnelling  would  be  far 
slower  for  than  for  because  of  the  greater  reduced  mass  (Bates,  Guest  and 

Kendall  1993). 

A  straightforward  (but  hitherto  overlooked)  argument  proves  that,  for  many 
polyatomic  ions,  the  crossings  of  potential  curves  that  had  been  presumed  to  be 
instrumental  in  the  dissociative  recombination  process  cannot  exist.  For  instance, 
in  the  dissociative  recombination  of  a  polyatomic  hydride  ion  leading  to  the  release 
of  an  H  atom,  a  single-electron  radiationless  transition  is  always  involved  and  tins 
precludes  a  crossing  of  potentials.  And  yet  for  every  such  case  on  which  experimental 
data  is  available  this  process  is  rapid.  Thus,  rapid  dissoc iative  recombination  without 
a  crossing  of  potentials  is  common.  The  experimental  results  may  be  accounted  for 
by  the  new  multi-step  cousin  of  the  original  mode  of  dissociative  recombination. 
This  cousin  is  more  common  in  nature  than  the  curve-crossing  form  of  dissociative 
recombination  (Bates  1993a,  1993d). 

A  considerable  fraction  of  any  molecular  product  of  the  dissociative  recombination 
of  a  polyatomic  ion  has  been  proven  to  be  normally  vibrationally  excited.  The  usual 
main  source  of  the  excitation  is  the  impulsive  force  between  the  products.  In  the  case 
of  a  monohydride,  the  vibrational  excitation  due  to  the  impulse  acting  on  its  heavy 


7 


atom  is  inappreciable,  but  that  due  to  the  much  weaker  impulse  acting  on  its  H  atom 
may  be  great.  Another  source  is  rooted  in  the  relic  of  the  normal  vibrational  modes 
of  the  ion.  This  causes  the  initial  vibrational  wavefunction  of  the  product  to  differ 
from  the  wavefunction  of  its  zero-th  vibrational  level,  even  in  instances  where  the 
corresponding  equilibrium  separations  force  constants  of  the  the  ion  and  the  product 
are  virtually  the  same  (as  they  commonly  are).  If  the  ion  is  bent,  as  for  example  in 
OCOH  +  ,  an  impulsive  couple  may  act  on  one  or  more  of  the  products  and  cause  high 
rotational  excitation  (Bates  1993b). 

3.  Tidal  termolecular  ionic  recombination 

Experimental  measurements  by  Mezyk  et  a 1  (1989,  1991)  have  shown  that  the  rate 
coefficient  for  the  process 

Xe  +  +  Cl-  +  Xe  - ♦  Xe  Cl*  +  2 X e  (2) 

is  enhanced  relative  to  that  given  by  the  universal  curve  that  describes  termolecular 
recombination  of  structureless  atomic  ions.  This  is  attributed  to  Cl-  exciting  internal 

-L 

modes  of  Xe2  during  passage  through  the  perihelion  part  of  the  orbit.  We  have 
simulated  the  electrostatic  tidal  effect  by  the  Monte  Carlo  technique  and  have  obtained 
excellent  agreement  with  the  experimental  data.  It  is  reasoned  that  the  tidal  action 
cannot  enhance  the  rate  of  termolecular  recombination  of  small  valence-bound  ions 
like  and  02  appreciably.  However,  its  influence  is  pervasive  because  complex 

ions  are  commonly  present  when  the  ambient  gas  density  is  enough  for  termolecular 
recombination  to  be  of  interest.  For  example,  is  the  main  positive  ion  if  the  ambient 
gas  is  pure  oxygen.  Enhanced  mutual  neutralization  may  obscure  the  experimental 
evidence,  but  may  be  disregarded  in  the  processes 

H30+  (H20)n  +  NO ^  +  {He,  Ar}  — >  neutrals.  (3) 

Experimental  research  by  Lee  and  Johnsen  (1989)  has  shown  that  the  rates  of  these  two 
processes  are  enhanced  relative  to  the  universal  curve.  The  only  possible  explanation 
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for  that  enhancement  seems  to  be  tidal  action.  A  paper  describing  our  preliminary 
computations  was  published  (Bates  and  Morgan.  1990),  and  a  paper  containing  a  more 
detailed  discussion  and  the  results  of  our  extensive  computations  has  been  published 
(Morgan  and  Bates,  1992b). 

4.  Radiative  and  termolecular  association 

Earlier  calculations  on  association  in  CH^"  +  H2  and  CD^  +  D2  collisions  have  been 
extended  and  refined  (Bates,  1992b).  It  was  concluded  that  there  are  serious 
discrepancies  between,  on  the  one  hand,  recent  measurements  made  by  Gerlich  and 
Kaefer  (1989)  with  an  80  K  low-pressure  ion-trap,  and,  on  the  other  hand,  theory  and 
measurements  made  by  Smith  et  a 1  (1982)  with  a  SIFT  apparatus  at  80  K,  0.4  torr 
and  measurements  made  by  Barlo\v  et  &l  (1984)  using  a  13  K  ultra-low-pressure  trap. 
No  explanation  could  be  offered. 

5.  Nightglow  emissions  from  oxygen  in  the  lower  thermosphere 

Several  independent  trains  of  reasoning  have  led  to  the  conclusion  that  02(c1E“) 
is  the  precursor  for  the  emission  of  the  green  line  of  atomic  oxygen  (Bates.  1992c). 
We  have  shown  that  the  measured  altitude  profile  of  the  volume  emission  rate  of 
the  Hertzberg  I  system  is  reproduced  by  using  the  expected  rate  coefficients  for 
the  production  of  02(A3D^)  by  termolecular  association  and  rate  coefficients  for 
quenching  by  atomic  and  molecular  oxygen  that  accord  with  laboratory  evidence. 
Contrary  to  the  common  belief,  we  have  inferred  that  termolecular  association  does  not 
contribute  directly  to  the  Atmospheric  system  because  the  02  (b  1  E  + )  it  generates  is 
mainly  in  high  vibrational  levels  and  deactivation  of  these  by  interchange  of  electronic 
energy  with  normal  oxygen  molecules  is  not  rapid  enough  to  prevent  deactivation  to 
02  (a  1  A  )  or  02(X3E”)  from  occurring.  We  have  shown  that  the  prime  source  of 
the  02  (b  1  Hg  u=0)  is  the  excitation  of  the  02(X3E~  u=0)  collision  partner  involved 
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in  02  (A  3  £*  )  ,  00(A/3Au)  and  02  (c  1  Eu  )  deactivation.  This  source  can  ac count  for 
the  measured  altitude  profile  of  the  volume  emission  rate  of  the  Atmospheric  system. 

6.  Nightglow  emissions  from  atomic  oxygen  in  the  ionosphere 

The  dissociative  recombination  quantum  yields  /('  S)  an  d  /(*  D)  that  h  ave  been 
deduced  by  Sobral  et  a 1  (1992)  from  the  nocturnal  rocket  measurements  on  the 
volume  emission  rates  V(5577)  and  1^(6300)  in  the  equatorial  F-region  have  been 
analysed  to  see  if  they  are  consistent  with  the  results  of  the  a b  initio  computations 
of  Guberman  and  Giusti-Suzor  (1991)  on  dissociative  recombination  (Bates.  1 9 9 2 d ) . 
In  the  case  of  the  557.7  nm  emission,  the  only  serious  problem  we  have  found  is  that 
the  rocket  f(lS)  values  are  much  larger  than  would  be  expected  from  the  ab  initio 
computations  in  the  region  below  about  220km.  The  most  interesting  feature  of  the 
apparent  /(*D)  altitude  profile  that  we  have  calculated  is  that  it  is  markedly  influenced 
by  the  transport  of  the  suprathermal  0(JD)  atoms  produced  by  the  dissociative 
recombination. 

7.  Nightglow:  the  source  of  the  green  continuum 

By  analysing  observations,  Barbier  et  a  1  (1951)  established  that  there  is  a  continuum 
in  the  green  region  of  the  nightglow  spectrum.  The  accepted  theory  (Krassovsky  1951) 
is  that  it  arises  from  the  N02  air  afterglow  from  the  composite  process 

NO  +  O  +  M  - >  N02  +  M  +  hi/  (4) 

with  the  NO  being  regenerated  by  the  fast  reaction  between  N02  and  0.  However, 
by  analysis  of  the  observational  evidence,  we  have  shown  that  this  standard  theory  is 
untenable.  Furthermore,  we  have  examined  and  eliminated  all  of  the  conceivable 
“stock”  sources  of  the  continuum  emission.  The  observational  evidence  points 
unambiguously  to  a  novel  source  of  the  continuum  emission,  namely,  metastable 
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oxygen  molecules  colliding  with  ambient  gas  molecules  and  thereby  forming  complexes 
that  dissociate  by  allowed  radiative  transitions. 

Let  O2  denote  an  oxygen  molecule  in  one  of  the  six  metastable  states  (a  1  A? ,  b  1  , 

c  1  ~ .  A'  3Au,  and  5ng)  that  are  populated  by  termolecular  association.  We 

propose  that  a  collision  between  O2  and  another  molecule  or  atom  is  intimate  enough 
to  form  a  complex  that  dissociates  radiatively  by  an  allowed  transition.  Inspection 
of  the  potential  curves  of  diatomic  molecules  shows  that  conditions  for  an  allowed 
transition  are,  in  many  instances,  created  when  a  metastable  atom  approaches  a 
normal  atom.  This  must  also  be  true  of  molecules,  with  the  vital  difference  that 
the  duration  of  the  intimate  collision  may  be  much  longer  because  of  the  high  energy 
density  of  the  rovibrational  levels  of  the  complex.  In  general  the  radiative  transition 
ensures  the  dissociation  of  the  complex,  so  that  the  emission  forms  a  continuum.  If, 
for  instance, 'the  collision  partner  is  a  normal  oxygen  molecule,  the  proposed  sequence 
is  then 

°2  +  02  (X  3S~)  O;  (5) 

04  ►  O2  +  O2  +  hi/ 

the  primes  being  reminders  that  the  products  are  not  necessarily  in  the  ground 
electronic  state. 

Using  rocket  data  (Greer  et  a 1  1986)  we  find  that  if  the  collision  partner  is 

02  (X  3n?  )  or  N2  (X  'S^")  the  emission  would  peak  near  96  km  and  that  the  thickness 
of  the  luminous  layer  (defined  as  the  vertical  distance  between  the  levels  where  the 
intensity  is  half  the  maximum)  would  be  10km.  These  values  are  in  satisfactory 
agreement  with  the  rocket  measurements  of  McDade  et  a 1  (1986). 

To  account  for  the  observed  nightglow  intensity,  the  maximum  volume  emission 
must  be  around  1  X  103  cm3s_1  in  the  visible  region  and  5  X  103  cm3s-1  in  the  infrared 
region.  We  have  estimated  that  such  values  can  be  attained  by  our  proposed  new 
mechanism,  although  laboratory  measurements  are  needed  to  check  the  estimated 
rates  that  we  have  used. 

This  work  has  been  published:  see  Bates  (1993c). 


11 


8.  The  ionosphere 


A  history  of  some  aspects  of  ionospheric  research  has  been  written  (Bates,  1993e). 
Seriously  mistaken  beliefs,  for  example  that  negative  ions  are  the  arriers  of 

negative  charge  in  the  E-layer,  have  been  widely  held.  lonograms  leu  perception 

on  the  diurnal  variation  of  iV,  the  number  density  of  free  electrons,  that  forced 
theorists  to  propose  that  electrons  and  molecular  positive  ions  disappear  by  a  process 
that  is  many  orders  of  magnitude  faster  than  had  earlier  been  inferred  for  any  other 
loss  process  envisaged.  The  proposal  regarding  recombination  is  correct,  1  never, 
incoherent  scatter  radar  studies  have  given  that  N  in  the  E-layer  has  an  asymmetry 
with  respect  to  noon  that  is  in  the  opposite  sense  to  that  deduced  from  ionograms. 
and  this  behaviour  has  been  satisfactorily  reproduced  by  numerical  simulation.  The 
harmony  is  not  complete,  good  reason  not  having  been  given  for  neglecting  the  results 
of  the  ionogram  research. 

9-  Heavy  carbon  dioxide  enhancement  in  the  stratosphere 

The  commonly  accepted  source  of  the  heavy  C02  enhancement  in  the  stratosphere  is 
transfer  from  the  heavy  03  enhancement.  It  is  reasoned  (Bates,  1991a)  that  another 
source  is  provided  by  the  effect  of  the  avoided  crossing  of  the  03(1  1  B2)  state  that  is 
reached  by  absorption  in  the  Hartley  band,  and  we  have  shown  that  the  sources  are 
equal  if  the  heavy  03  enhancement  is  by  about  20  per  cent. 

10.  Heavy  ozone  in  the  stratosphere 

Flips  in  which  one  of  the  end  atoms  in  an  energized  complex  like  160  160  180  *  ceases 
to  be  effectively  bound  to  the  central  atom  and  becomes  effectively  bound  to  the  other 
end  atom  (so  that  0  is  transformed  into  the  central  atom)  have  been  treated.  Using 
the  Monte  Carlo  method  to  follow  the  motion  of  three  oxygen  atoms  on  the  ground- 
state  potential-energy  surface  of  ozone,  it  has  been  shown  that  such  flips  do  indeed 
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occur  and  that  their  frequency  relative  to  the  dissociation  frequency  is  consistent  with 
the  available  experimental  results  (Morgan  and  Bates  1992a). 


11.  Photoionization  of  atomic  nitrogen 

A  study  of  the  photoionization  of  atomic  nitrogen  was  made  (Bell  et  a l  1992)  by 
extending  the  work  of  Bell  and  Berrington  (1991)  to  a  consideration  of  photoionization 
of  the  six  excited  states  2s22p3  2D°,  2s22p32P°,  2s2 2p2 (3 P)3s  4 P ,  2s22p2(3P)3s  2P. 
2s2p4  4P  and  2s2  2p2  (3  P)3p  2  S° .  The  i?-matrix  method  was  used  and  thirteen 
target  states  of  N  +  have  been  included  in  the  calculation:  2s22p2  3P,  *D,  'S; 
2s2p3  5S°,  3D°,  3P°,  1D°;  2s22p3s  1  P°,  3P°;  2s2p33S°;  2s22p3p  1  P,  3D ;  ‘2s2p3  1P°. 
Each  of  these  thirteen  N+  states  was  represented  by  configuration  interaction  wave- 
functions. 

The  present  results  for  the  2s22p3  2D°  and  2s22p3  2P°  states  can  be  compared 
with  the  earlier  less  sophisticated  work  of  Zeippen  et  a 1  (1980).  The  close  agreement 
between  the  two  calculations  is  significant  in  that  it  demonstrates  that  the  R-matrix 
calculations  have  converged  and  that  the  cross  section  is  now  known  to  an  accuracy 
of  better  than  a  few  per  cent.  The  agreement  is  particularly  satisfactory  in  the  region 
of  the  2s2p4  2P  resonance  and  the  background  cross  section  in  general  in  both  cases. 

The  cross  sections  for  2s2  2p2  (3  P)3s  4  P  and  2s22p2  (3  P)3s  2  P  are  notable  for  an 
increase  in  the  value  of  the  cross  section  at  the  3s 3  P°  threshold  —  the  increase 
corresponding  to  the  ionization  of  the  2p  electron  into  the  continuum  with  the  3s 
electron  behaving  like  a  ‘spectator’  electron.  A  similar  increase  is  found  for  the 
2s2p4  4P  case.  No  experimental  data  exists  for  any  of  these  cases. 

For  the  2s22p2(3  P)3p  2S°  state  the  ‘kink’  or  shallow  minimum  in  the  cross  section 
(associated  with  a  Cooper  minimum)  is  confirmed.  However,  considerable  divergence 
of  the  present  work  for  the  asymmetry  parameter  for  this  state  occurs  in  comparison 
with  earlier  calculations  (Theodosiou  1988,  Nahar  and  Mason  1989,  and  Mason  1988). 
Although  the  single  experimental  point  (Pratt  et  a 1  1987)  is  insufficient  to  support 


13 


any  theoretical  result,  we  believe  that  the  present  calculation  is  the  most  accurate 
presently  available. 


12.  Photoionization  of  atomic  oxygen 


The  ^-matrix  method  has  been  used  to  calculate  the  cross  sections  for  photoionization 
of  ground-state  atomic  oxygen,  leading  to  the  production  of  several  possible  final-ion 
states: 


O  (ls22s22p4)3P  +  hv 


0+  (ls22s22p3)4S°,  2D°,  2P°  +  e“ 
0+  (ls22s2p4 )  4  P.  2P  +  e“. 


(6) 


Data  on  these  processes  are  of  fundamental  importance  in  the  quantitative  under¬ 
standing  of  the  ionosphere. 

Details  of  the  i?-matrix  codes  for  photoionization  calculations  have  been  given 
by  Berrington  et  a 1  (1987),  with  the  associated  theory  described  by  Seaton  (1987). 
The  present  calculation  is  done  in  XS-coupling  and  negiects  relativistic  effects.  The 
states  of  0+  included  in  the  wavefunction  expansions  are:  2s22p3  4S°,  2D°,  2P°, 
2s2p4  4P,  2D,  2P;  2s22p23s4P,  2P,  2D,  2S.  Each  of  these  ten  final-ion  states  was 
represented  by  configuration-interaction  type  wavefunctions  obtained  by  employing 
the  Is,  2s  and  2p  orbitals  of  Clementi  and  Roetti  (1974)  and  the  3s,  3p  and  3d  orbitals 
of  Bell  et  al  (1989),  The  R- matrix  radius  was  taken  to  be  12.4  aQ,  and  20  continuum 

orbitals  were  used  for  each  of  the  required  angular  momenta  (0 . 5).  The  set  of 

correlation  functions  was  comprised  of  all  configurations  arising  from  2s2p3.  2s22p2 
and  2p4  plus  two  electrons  from  the  given  set  of  bound  orbitals.  The  cross  sections 
were  calculated  in  the  length  formulation,  and,  in  the  light  of  previous  work  (Bell 
et  ai  1987),  we  believe  that  the  our  total  photoionization  cross  sections  should  be 
accurate  to  a  few  percent.  We  thus  claim  that  the  present  work  on  photoionization 
of  ground-state  atomic  oxygen  is  the  most  reliable  calculation  yet  performed.  It  has 
been  published:  see  Bell  and  Stafford  (1992). 

Our  results  for  photoionization  in  which  a  2p  electron  is  removed,  producing 
0  +  (2s22p3  4S°),  0+(2s22p3  2D°)  or  04’(2s22p3  2P°),  are  not  too  disparate  from  other 
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calculations.  Thus,  for  photoionization  that  yields  0  +  (2s22p3  4S°)  our  results  are  in 
close  agreement  with  those  of  Pradhan  (1978)  and  Kirby  et  al  (1979),  except  that  for 
wavelengths  below  360  A  our  data  decrease  more  rapidly  than  those  of  Kirby  et  al  . 
Similarly,  our  results  for  photoionization  that  produces  0+(2s22p3  2D°)  agree  closely 
with  the  data  of  Kirby  et  al  (1979),  although  they  are  significantly  above  the  values 
reported  by  Pradhan  (1980).  The  agreement  between  our  work  and  other  work  is 
not  so  good  for  photoionization  that  produces  0  (2s  2p  P  ).  Our  results  are  below 

those  of  Kirby  et  al  (1979)  and  of  Pradhan  (1980). 

For  photoionization  in  which  a  2s  electron  is  removed,  producing  0+(2s2p4  4P) 
or  0+(2s2p42P),  our  results  differ  dramatically  from  previously  reported  data.  The 
present  results  have  led,  for  the  first  time,  to  a  theoretical  fit  in  agreement  with  the 
rocket  data  on  the  834  A  airglow  (Link  et  al  1993). 


13.  Electron-molecule  collisions 

Extensive  calculations,  using  the  R- matrix  method,  have  been  made  on  electron 
scattering  by  oxygen  molecules.  The  /2-matrix  method,  as  applied  to  electron- 
molecule  scattering,  was  first  described  by  Burke  et  al  (1977),  with  the  essential 
techniques  set  out  by  Gillan  et  al  (1987).  Electron  excitation  of  02  to  the  and 

61E^'  states  has  been  studied  by  the  /2-matrix  method  by  Noble  and  Burke  (1986), 
but  they  considered  only  the  three  lowest  electronic  states  of  the  the  target. 

In  the  first  calculations  of  the  present  study  (Noble  and  Burke  1992),  the  lowest 
nine  states  of  02  were  included  within  the  /2-matrix  framework.  Electronic  excitation 
cross  sections  for  transitions  to  the  alAg,  and  ‘6eV’  states  were  calculated 

for  impact  energies  up  to  15  eV .  Important  contributions  from  a  IIU  resonance,  in 
addition  to  the  well-known  resonance,  were  found  to  modify  significantly  previous 
theoretical  predictions.  Also,  integral  elastic  cross  sections  were  obtained  for  the  first 
time  for  a  complicated  open-shell  system  using  an  ab  initio  method . 

In  the  second  phase  of  the  study  (Higgins  et  al  1994),  the  earlier  work  by 
Noble  and  Burke  (1992)  has  been  extended.  The  same  nine  states  of  02  have 
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been  included,  but  the  Cl  wavefunctions  describing  thou,  states  were  improved  by 
including  additional  configurations.  Also,  a  broad  range  of  values  of  R.  the  inter- 
nuclear  separation,  have  been  considered  —  the  earlier  work  considered  the  equilibrium 
value  only,  R  =  Re  %  2.3  aQ.  The  cross  sections  have  been  calculated  for  elastic 
scattering  and  for  electronic  excitation  from  the  X  3E~  ground  state  to  the  a1.^, 
b1!^  and  the  ‘6eV’  states  and  compared  with  experimental  data  whenever  possible. 
The  positions  and  widths  of  the  resonances  in  the  II  ^  and  IIU  scattering  symmetries 
have  been  calculated  as  functions  of  R  for  1.85  aQ  <  R  <  2.90  a0,  and  their  effect  on 
the  cross  sections  determined.  The  results  show  the  importance  of  using  a  good  target 
representation  in  order  to  obtain  consistent  results  as  functions  of  R.  The  value  of  R 
is  also  found  to  have  a  profound  effect  on  the  scattering  cross  sections,  particularly  in 
the  resonant  2IIS  and  2IIU  symmetries. 

The  data  produced  by  the  iZ-matrix  calculation  of  Noble  and  Burke  (1992)  has  been 
used  to  study  the  change  in  polarization  of  transversely  polarized  electrons  scattered 
by  randomly-oriented  oxygen  molecules  in  their  ground  X3£~  state.  The  theory  of 
how  this  information  can  be  calculated  from  the  T-matrix  elements  produced  in  the 
iZ-matrix  calculation  has  been  described  by  Bartschat  and  Madison  (1988).  Changes  in 
polarization  can  be  due  to  exchange  interactions  or  spin-orbit  interactions,  the  former 
being  expected  to  be  more  significant  for  a  light  molecule  such  as  02.  Hence,  our 
study  should  cast  light  on  the  important  issue  of  the  relative  strength  of  the  exchange 
interaction  in  low-energy  electron-molecule  collisions. 

We  have  found  that  the  polarization  fraction  for  elastic  collisions  from  02  (X  3E ~) 
is  close  to  unity  for  all  angles  and  energies  between  1  eV  to  15  eV.  This  confirms  the 
experimental  findings  of  Hegemann  et  a 1  (1990),  and  the  results  from  the  earlier, 

simpler  calculation  of  da  Paixao  et  a 1  (1992).  The  2ng  resonance  produces  sharp 

structure  for  energies  below  1  eV .  Structure  due  to  the  2I1U  resonance  at  8eV  can 
be  seen  at  small  and  large  scattering  angles,  but  the  effect  is  small  and  probably  not 
measurable.  For  energies  between  10  eV  and  15  eV,  our  results  may  be  sensitive  to  the 
resonance,  and  agreement  with  experiment  in  this  energy-range  may  be  improved 
by  more  extensive  scattering  calculations.  Our  results  for  the  polarization  change  in 
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collisions  that  excite  the  oxygen  molecule  are  the  first  to  be  reported.  We  have  found 

a  wide  variation  in  the  polarization  fraction  as  a  function  of  angle  and  as  a  function 

<2 

of  energy,  a  manifestation  of  the  strong  influence  of  the  Ilu  resonance.  This  work  has 
been  published:  see  Fullerton  et  al  (2994). 

14.  Electron-impact  ionization  of  metastable  inert  gases 

Total  and  single  differential  cross  sections  for  electron  impact  ionization  of  helium 
atoms  have  been  calculated  using  the  R-matrix  method  of  Bartschat  and  Burke  (1987). 
In  addition  to  the  ground  state,  ionization  from  the  two  metastable  states  has  been 
investigated,  and  the  possibility  of  producing  an  excited  helium  ion  has  been  included: 

He  ( Is2  ) 1 S  'j 

He  ( ls2s)  1 S  >  - ►  e(klll)  +  He+(nf)  +  e(k2l2)  •  (7) 

He  (ls2s)  3S  J 

where  kt  and  denote  the  wave-numbers  and  angular  momenta  when  the  electrons 
are  described  by  partial-wave  expansions.  Six  final  states  of  He+  were  included:  Is, 
2s,  2p,  3s,  3p  and  3d.  The  2^-pole  components,  with  A  <  4,  of  the  ef  +  He  interaction 
were  included,  giving  five  possible  He-continuum  symmetries.  Due  to  the  inclusion  of 
states  of  the  final  ion  with  n  =  3  and  to  the  initial  state  being  excited,  it  was  necessary 
to  use  a  large  value  of  the  R-matrix  boundary,  necessitating  a  large  number  (25)  of 
continuum  orbitals.  The  calculated  energies  of  the  initial  states  of  He(l  1 S ) ,  He  (2  lS) 
and  He(23S)  indicates  that  the  bound  states  generated  by  the  R-matrix  method  are 
reliable. 

The  results  have  been  compared  with  experiments  and  other  theoretical  data.  The 
discrepancy  between  our  results  and  experiment  for  total  cross  sections  suggests  that 
there  is  a  need  to  include  the  effects  of  exchange  between  the  “fast”  election  and 
the  electrons  of  the  target.  The  discrepancies  between  our  results  and  those  from 
other  theoretical  approaches,  especially  the  close  coupling  results  of  Rudge  (1988), 
demonstrate  the  sensitivity  of  the  calculated  cross  sections  to  the  wavefunctions  of 
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the  initial  bound  state  and  the  final  continuum  states,  and  perhaps  also  to  the  use 
of  distorted  waves  rather  than  plane  waves  for  the  "fast'’  electron.  We  take  this 
as  vindicating  our  use  of  the  R-matrix  formulation,  even  though  the  excitation- 
autoionization  process,  for  which  the  R-matrix  formulation  is  particularly  adept,  is 
not  significant  in  the  present  study.  The  close  coupling  in  the  final  continuum  states 
is  crucial  for  excitation  of  final-ion  states,  such  as  He+(2p),  that  cannot  be  produced 
by  simple  removal  of  one  electron. 

This  work  on  helium  has  been  published:  see  Raeker  et  al  (1994).  It  has  provided 
a  valuable  test  of  the  method  for  excited  targets,  enabling  us  to  proceed  to  consider 
ionization  of  metastable  Ne  and  Ar. 
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Final  Remark 


Sir  David  Bates  died  suddenly  on  January  5,  1994.  The 
news  of  his  death  was  received  with  great  sadness  and  sense 
of  loss  by  his  colleagues  at  The  Queen’s  University  of  Belfast 
and  by  his  many  friends  and  collaborators  throughout  the 
international  science  community.  This  report  has  been 
completed  by  Sir  David’s  colleagues  ~.t  Belfast,  who  feel  sure 
that  he  would  have  wanted  them  to  express  his  gratitude 
to  the  Air  Force  Office  of  Scientific  Research  for  their 
support,  which  has  extended  over  many  years.  This  support 
enabled  Sir  David  to  pursue  his  research  with  undiminished 
enthusiasm  and  vigour  to  within  days  of  his  death. 
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